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Abstract 

We present results for the reference scale r$ in SU(3) Lattice Gauge Theory for /3 = 
6/<?o i n t ne range 5.7 < (3 < 6.57. The high relative accuracy of 0.3-0.6% in r^/a 
was achieved through good statistics, the application of a multi-hit procedure and a 
variational approach in the computation of Wilson loops. A precise definition of the 
force used to extract ro has been used throughout the calculation which guarantees 
that ro/a is a smooth function of the bare coupling and that subsequent continuum 
extrapolations are possible. The results are applied to the continuum extrapolations 
of the energy gap A in the static quark potential and the scale L max Ao used in the 
calculation of the running coupling constant. 
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1 Introduction 



Monte Carlo calculations of lattice QCD in the quenched approximation and in the 
SU(3) Yang-Mills theory have reached considerable precision. For example, the accuracy 
of hadron mass calculations using the Wilson action is quoted to be well below 1% at 
finite values of the lattice spacing Q , and the running of the coupling in a specific non- 
perturbative scheme is known to a precision of around 1% over energy scales varying by 
two orders of mag nitude In the latter case, the continuum limit was taken. 

It was already noted in Q that in comparison to such a level of accuracy, a low- 
energy reference scale in the Yang-Mills theory is known with much worse precision, 
despite the fact that such a scale is very important for the analysis of the results. In 
the second example of the running coupling, a low-energy reference scale is required in 
order to determine the overall momentum scale for which the coupling was computed 



(see Sect. 3.2). In the first example, a gluonic reference scale is very useful in order to 
assess the size of lattice artifacts present in different discretizations of QCD || and to 
provide a check on the continuum extrapolations presented in Jl], which have been per- 
formed using hadronic scales. The reason is that, in the quenched approximation, the 
leading cutoff effects in gluonic observables are proportional to a 2 , only. For instance, 
for Wilson fermions without improvement, the ratio of a hadronic and a gluonic mass 
scale approaches the continuum limit with an asymptotic rate proportional to the lat- 
tice spacing a, where the linear a-dependence originates purely from the hadron mass. 
Obviously, good precision in the gluonic reference scale is desirable for such applications. 

In a reference scale, ro, was introduced for such purposes. This length scale is 
defined in terms of the force, F(r), between external static charges in the fundamental 
representation. It is the solution of 

r 2 F(r ) = 1.65. (1.1) 

The constant on the right hand side was chosen such that ro has a value of approximately 
0.5 fm in QCD [||. In comparison to glueball masses and the string tension, it is much 
easier to compute with controlled errors [J 

However, for gauge group SU(3) and the Wilson plaquette action, ro has been 
computed only with modest precision so far. A recent new effort concentrated on 
the region of small (3 / large lattice spacing. In this paper we compute ro down to 
small lattice spacings a ~ 0.04 fm, covering the whole range 0.17 fm < a ^ 0.04 fm and 
improving the precision by a factor up to five. This gain in precision is only to a 
small part due to the use of a modern parallel computer (APE100), but results mainly 
from the application of the known methods of variance reduction Q and a variational 
calculation ||||], as employed already in ||] for gauge group SU(2). 

x The reason is that glueball correlation functions are quite noisy; the string tension refers to the 
asymptotic force at infinite distance, which is difficult to extract without additional assumptions about 
the leading corrections to the large-distance behaviour. 
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Our results are given in the form of a table and also an interpolating fit function 
which describes ro with a precision between 0.3% (at j3 sa 5.7) to 0.6% (at « 6.57). 
At the level of such precision, one can expect that the results depend on the chosen 
discretization of the force F(r). One and the same definition has to be used at all 
values of (3 in order to guarantee a smooth approach to the continuum limit. The 



procedure applied in j|] and [10| does not correspond to a precise definition of the force 
satisfying this criterion. Indeed, our results differ slightly from the ones reported in 
@,|0|. The differences are only 2.3 ± 0.9% at = 5.7, 0.9 ± 0.5% at = 5.85 and 



— 1.2 ± 0.6% at = 6.2. From our study of finite size effects we can exclude that these 
differences are due to different lattice sizes used in the various calculations. 

We apply our results for ro to the calculation and continuum extrapolation of the 
ratio L max /ro, which determines the scale in the running coupling and A-parameter 
calculations of ||. We also compute the gap between the ground state potential and its 
first excitation in the continuum limit. 

This paper is organized as follows. In the following section we describe the com- 
putation of ro in lattice units, starting with a recapitulation of the precise definition. 



Readers uninterested in these details will find the results in Sect. 2.4 and the applications 



to two continuum extrapolations in Sect. |3|. We finish with a brief conclusion. 

2 Determination of r 
2.1 Definition of the force 

The definition eq. (|bl|) is unique in the continuum. For finite lattice spacing, we have to 
specify which discretization of the force is to be used. Lattice artifacts - and therefore 
the precise values of r^/a - do depend on these details. However, it has become cus- 
tomary to replace the precise definition Q by a fitting procedure, with fits performed 
to potential values in the neighbourhood of ro [11,1^,|]. The motivation for such a 
procedure is to increase the statistical precision, but then it is not guaranteed that r^/a 
is a smooth function of 0, allowing for systematic continuum extrapolations. We will 
therefore return to a precise definition and demonstrate that the achieved statistical 
precision is in fact very good. 

We start from potential values, V(r), given along one fixed orientation d/|d| on 
the lattice. Here, d is a (three-dimensional) vector on the lattice, and r is an integer 
multiple of d. For example we may have d = (a, a, 0). The force is then defined as 

FiM = y %7 r - d > (2.1) 
1 _ G(r)-G(r-_d) 



47rrj |d 



_l /" d 3 A; U 3 j=i cos ( r jkj/a) 
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The particular choice of r\ ensures that the force is given by i*d(n) = 3ff 2 /(4vrrj )+0(g 4 ). 
To lowest order of perturbation theory, it eliminates the lattice artifacts exactly; they 
remain (probably quantitatively reduced) only in the higher 0(g 4 ) terms. The force 
Fd{ r i) as defined above is therefore called a tree- level improved observable. 

In the following, we choose d = (a, 0, 0) dropping this index on the force. This 
choice has also been made in and our data may therefore be directly compared with 
the results from that reference. Solving eq. (1.1) requires furthermore an interpolation 



of F(r\), which is easily done with small systematic errors ||. 

If the procedure to compute the scale ro, as outlined in this subsection, is applied 
consistently at all values of (5 considered, then r^ja will be a smooth function of the 
bare coupling. 

2.2 Wilson loop correlation matrix 

We now turn to the details of the computation of V(r). It is well known that an effective 
calculation of the potential starts from smeared Wilson loops. Our smearing operator 
S acts on the spatial components of the gauge fields via p 



SU(x,k) = V{U(x,k) + a J2 [U(x,j)U(x + aj,k)U\x + ak,j) + 

U\x - aj,j)U(x - aj, k)U(x + ak - aj,j)]} , (2.4) 

where V denotes a projection back into the group SU(3) and the SU(3) gauge field 
residing on a link from x to x + aft is denoted by U(x, /x). For different smearing levels 
I = 0, . . . , M — 1 we then construct smeared spatial links according to 

Ui{x,k) = S ni U(x,k). (2.5) 

For the time-like links in the Wilson loops, we apply the multi-hit procedure 0] in order 
to reduce the variance. These multi-hit averaged linksf] are denoted by U . At fixed r, 
an M x M correlation matrix of Wilson loops is then formed as 

Q m (t) = ( tr {^(0, ri)F(rl, A + t6)Vl{t6, rl + t6)V\o, i6)}} = C ml {t) , (2.6) 

where 

Vi(x,x + rl) = Ui(x,l)U l (x + al,l)...U l (x + (r-a)i,l), (2.7) 

V(x,x + t6) = U(x,0)U(x + a6,0)...U(x + (t-a)6,0). (2.8) 

These correlation matrices were computed on lattices generated by a hybrid over- 
relaxation algorithm [14,l]| with a mixture of N or ~ 1.5 ro /a over-relaxation steps per 



heat bath step. In each simulation, a measurement of Wilson loops is performed every 



2 We applied an average using 10 Cabibbo-Marinari iterations each of which consisted of updates in 
three SU(2) subgroups embedded in SU(3). 
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L/a 


P 


-^meas 


M 


L/a 





-^rncas 


M 


10 


5.7 


1000 


4 


24 


6.2 


400 


4 


12 


5.8 


3200 


4 


32 


6.4 


200 


3 


16 


5.95 


600 


4 


40 


6.57 


572 


1 


20 


6.07 


780 


4 











Table 1: Simulation parameters. iVmeas denotes the number of different gauge field configurations for 
which the correlation matrices were computed. For every configuration averages were computed over 
all (L/a) 4 points of the lattice and three permutations of the lattice axes. M is the number of smearing 
levels. 



20 x (iVor + 1) sweeps. The usual binning-jackknife procedure was used to compute the 
errors, and to control autocorrelations in the measurements. We further chose 

« = l n^U^) 2 . (2-9) 



2 ' ' 2 V a . 

The value we used for rii corresponds roughly to what was estimated to be the optimal 
smearing in j|. Other parameters are listed in Table [l]. Note that for [5 = 6.57 memory 
limitations restricted us to one smearing level (M = 1), and deviating from eq. 
we used no = 160 in this case. 



2.3 Ground state potential and gap 

Following [HH, the correlation matrices were analyzed using a variational method: take 
r and to fixed (to = a in practice) and solve the generalized eigenvalue problem, 

C(t)v a (t) = X a (t)C(t )v a (t) , (2.10) 

for the eigenvalues X a (t). These are related to the different potential levels by |J 

aV a = ln{X a (t)/X a (t + a)) + O^"^-^') . (2.11) 

Here, the potential levels V a are the eigenvalues of the (lattice) Hamilton operator in the 
sector of the Hilbert space that has the quantum numbers of the usual central potential 
V(r) = Vq. Besides the potential, we are also interested in the gap A = V\ — Vb, i.e. 
the energy difference between the ground state and the first excited state, since this gap 
controls the finite time corrections in the extraction of V(r) (see below). 
From eq. ( 2.11 ) we infer that A can be estimated by 

aA = ln(Ai(t)/Ai(t + a)) - ln(X (t) / X (t + a)) + CKe-^"^*) . (2.12) 

Checking our data for convergence at large t, we observed that the systematic corrections 
in eq. ( 2.12j ) are smaller than the statistical errors for t > |ro- This value of t was then 
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Figure 1: Estimates ro/a for various minimum time separations, t m in, in the fits eq. (2.14). (i increases 
from bottom to top and the values for different f3 have been shifted relative to each other for clarity. 
The ^-coordinate corresponds to the slowest varying finite-i correction term (as a function of i m m). For 
the gap A the value A « 3.3/ro has been used. Dashed error bands denote our final estimates. 



used for our estimates of A. We will return to these results in Sect. 3.1. Here we only 
note that A| r=ro ~ 3.3/ro in the continuum limit. 

As a first step in our analysis, the corrections to the ground state potential V(r) 



itself were studied using eq. ( 2.11 ) and found to be much smaller than in the case of the 
gap A. However, we were not satisfied with the statistical errors of ln(Ao(i)/Ao(i + a)). 
We therefore considered an alternative way of extracting V(r) by forming the projected 
Wilson loop correlation 

W(t) = v^C(t)v , (2.13) 

where vq is the eigenvector computed for t = to + a> in eq. ( |2.10|) . The projected 
correlation W{t) then turns out to be strongly dominated by the ground state. Next, 
W{t) was fitted to a single exponential 

W{t) w Ae- v W (2.14) 

for t ranging from t m \ n to t m3iX , usually setting £ max = i m i n + 3a. The parameters in these 
fits converge to the true values with corrections of order exp(— t m ; n A). To demonstrate 
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that these corrections are under control, we turn directly to the desired quantity rg. 



Using V(r) obtained from fits with fixed i m i n , we computed ro as outlined in Sect. 2.1 



We show in Fig. [j] the values of ro as a function of exp(— t m ; n A). Convergence is 
observed for large values of tmin- As our final values we took the fits with t m j n such that 
exp(— t m i n A) < 0.3. Let us summarize some important points about this analysis. 

• Like in the calculation of any energy value from Euclidean correlation functions at 
finite t, it is important to have an estimate for the gap in the channel considered. 
This is because time separations of the order of the inverse gap have to be reached, 
in order to judge whether correction terms are significant; 

• Due to the use of the variationally determined "wavefunction" vq, the corrections 
are most likely dominated by higher states and decay faster than exp(— tA) for 
the accessible values of t. Some evidence for this is seen in Fig. [l], in particular at 
larger values of f3; 

• In most cases the fits are statistically acceptable already for smaller t m [ n than that 
used for our final estimates; 

• Our estimates obtained directly from the application of eq. (2.11) for large t (start- 
ing from t = 3a) agree with our final estimates shown above; 

• The variance reduction Q] is essential to obtain estimates with moderate errors 
at large t. 

These observations make us confident that our estimates for the potential and ro rep- 
resent solid results. 



2.4 The scale r^/a 

The calculation of the force F(n) from the ground state potential and the subsequent 



interpolation of the force to extract r^/a follows the procedure outlined in Sect. 2.1 



2.4-1 Results; finite size effects 

The final values of r^/a are listed in Table ^. We note that the parameters (3 and 
L/a have been chosen such that L/r$ ~ 3.3, and hence the physical volume is kept 
constant over the whole range of (3 we considered. In order to check whether ro still 
has a significant dependence on L for L > 3.3ro, we carried out a series of calculations 
at (3 = 5.95 fixed but for several different values of L/a. The force F(r) for r in the 
vicinity of ro is shown in Fig. ||. One observes that, within errors, F(r) is independent 
of L for L/a > 16, which corresponds to L/r$ > 3.3. Therefore, we may take the values 
of ro listed in Table ^ as estimates in infinite volume. Table |2| also contains results from 
the literature at the same values of (3. One observes the small but significant differences 
quoted in the introduction. 
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(3 r /a L/r 



5.7 2.922( 9) 

5.8 3.673( 5) 
5.95 4.898(12) 
6.07 6.033(17) 
6.2 7.380(26) 
6.4 9.74 (5) 
6.57 12.38 (7) 




r /a 
2.990(24)" 



7.29(4)" 



ref. 



ig| 



Table 2: Results for ro/a. Our values and the respective volumes on which they have been computed 
are shown in columns 2 and 3. In the last four columns we show a comparison to other results found in 
the literature at the same /3 values. An asterisk indicates that the precise definition of ro has not been 
used. 




Figure 2: Finite size dependence of the force F(r). 



2-4-2 Parametrization 

We now describe the parametrization of our results for ro/a in Table [2] in terms of a 
smooth function of (3. This is meant to provide an interpolating formula, so that ro/a 
can be obtained at arbitrary values of (3 in the interval 5.7 < j3 < 6.57. 

A convenient starting point for the parametrization is the solution of the renormal- 
ization group equation for the bare coupling 



— = ^e-V^OSo 2 ) (6 o5 2)-6i/(2feg) e -ci9 2 -c 2 ^- 



(2.15) 



where b = ll/(4vr) 2 and b x = 102/(4tt) 4 are the universal one- and two-loop coefficients 
in the perturbation expansion of the (3 function, and the constant A is related to the 
A-parameter. The contributions containing c±, c%, . . . arise from higher order terms in 
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Figure 3: The data for ln(a/ro) (circles) and their representation as a polynomial in (3 (solid line). 

the perturbative f3 function. From the leading behaviour we infer 

(a/r )«e- /3 /( m °), (3 = 6/g 2 , (2.16) 

so that we attempt a phenomenological representation of ln(a/ro) in terms of a polyno- 
mial in (3. In order to avoid large cancellations among the fit parameters, it is advanta- 
geous to shift the value of /3, so that the intercept of the fit formula is contained in the 
interval of (3 values we considered. We choose the following ansatz 

ln(a/r ) = £ a fc (/? - 6) fe . (2.17) 

k=0 

A good description of our data is obtained for p = 3: 

In(o/ro) = -1.6805 - 1.7139 (J3 - 6) + 0.8155 {(3 - 6) 2 - 0.6667 (/? - 6) 3 , (2.18) 

and a comparison of this expression with the data points is shown in Fig. [31. The 
deviation between the curve and the values of ro/a is smaller than the statistical ac- 



curacy of the data. We thus take eq. ( 2.18; ) as our representation of ro/a in the range 



5.7 < (5 < 6.57. When ro/a is evaluated using eq. ( 2.18| ), a relative uncertainty of 0.3% 



at (3 = 5.7 should be assigned to the result, growing linearly to 0.6% at /3 = 6.57. This 
level of precision roughly corresponds to the statistical uncertainty in the data points. 
All results found in the literature, which use precisely our definition of tq [17,2,12.16] 
are described by eq. ( |2.18 ) within their statistical accuracy. 
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Figure 4: Continuum extrapolation of ro A. 



Another functional form of the parametrization is given directly by eq. ( |2.15| ). 



Since the three- loop contribution c\ has been calculated in [18], the fit parameters are 
A, C2, C3, .... In order to represent the data at the level of precision quoted above, four 
fit parameters are needed as in the phenomenological fit eq. ( |2.18| ). The fitted values of 
the parameters Cj turn out to be too large to allow for their interpretation as (effective) 
perturbative coefficients - as was expected from the well-known failure of "asymptotic 
scaling" in the bare coupling. We conclude that eq. ( 2.15| ) does not lead to a superior 
representation of the data. 



3 Continuum extrapolations 

In this section we want to give two examples of continuum extrapolatations with the 
scale set by ro. The first one is the gap A. To our knowledge this is the first time the 
location of an excited state is computed in the continuum limit of the 4-dimensional 
Yang-Mills theory. 



3.1 The potential gap 

The gap A was computed as a function of r as described in the previous section. Reliable 
estimates could be obtained for (3 < 6.2 only, since for the larger lattices the basis of 
operators (i.e. the parameter M) was too small. For illustration, we concentrate here 
on A| r=ro , noting that for smaller r the gap increases, while for much larger distances 
it becomes difficult to determine. To obtain the gap at r = ro, an interpolation in r 
must be performed, which can be done easily since the dependence on r is rather weak. 
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We then computed the dimensionless combination r^A. The errors of this quantity 
are much larger than for ro alone, namely around 2-5%. Nevertheless, as shown in 
Fig. |], lattice artifacts are clearly visible. They amount to more than 20% at f3 = 5.8 
(a/ro = 0.27) and decrease to around 12% at (3 = 5.95 (a/ro = 0.20). 

A continuum extrapolation can be performed using a linear fit in the leading cor- 
rection term (a/ro) 2 . Our graph shows the result obtained after discarding the data 
point at the largest value of the lattice spacing. The extrapolation yields 

r A| r=ro = 3.3(1) (3.1) 

in the continuum limit, and a similar number (with larger error) is obtained if one 
removes another data point from the fit. 



We note that bosonic string models (see e.g. |19,2C,2l[|) predict A = 2ir/r for large 
distances r. However, the fact that eq. ( |3.1| ) is almost a factor two smaller than this 
universal result should not be taken as evidence that an effective bosonic string model 
does not describe the QCD-string. Rather we have to remember that the universality 
of string models holds only for large r. For instance, for the Nambu-Goto action, the 
order 1/r 3 terms of the excitations are numerically very large at r = ro ||22|| . In contrast, 
the universal — 7r/(12r) term receives smaller order 1/r 3 corrections in that particular 
model, and this term has been observed for SU(2) Yang- Mills theory ||. We will return 
to an analysis of the 1/r-term in the potential for gauge group SU(3) in the future. 
Similarly to what we find, ref. [21] observes that a string model describes the Wilson 
loops of the three-dimensional Z2 gauge theory only for large loop sizes. 



3.2 The ratio L max /ro 

We now turn to discuss how to set the energy scale L max in the computation of the 
running coupling. For motivation let us recall the significance and the definition of L max 
p|, p3| , p^ ,^| . An important point in the quantitative and even qualitative understanding 
of QCD is to connect low-energy observables like hadron masses to the running coupling 
measured at high energies. If one adopts a suitable intermediate renormalization scheme, 
one can first connect the perturbative high energy region non-perturbatively to the 
coupling at some low-energy scale, q = 1/L max . In fact, as a result of being able to start 
deep in the perturbative region, one can compute the QCD A-parameter in units of L max 
||. Since the continuum limit can be taken in this calculation the result is universal, 
i.e. independent of the lattice action. A second step then involves the relation of L max 
to a low-energy scale like ro- Our aim is to perform this second step in this section. 



3.2.1 Schro dinger functional scheme 

A possible intermediate renormalization scheme is the Schrodinger functional scheme 
p5|| . The computation of the product L max Aj^g for quenched QCD will be described in 
detail in [26], and a summary has already been given in ||. Here we compute L max /ro, 
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which will ultimately allow to express A^jg in physical units. Note that quenched QCD is 
completely equivalent to the pure gauge theory in this context, since we are considering 
observables that do not involve fermion fields. 

In the Schrodinger functional scheme one considers QCD with specific Dirichlet 
boundary conditions in time (2f| (at xq = and xq = L), and the running coupling is 
defined in this scheme through an infinitesimal variation of the boundary values. For 
most details we refer to [Q], but for the following we need to briefly discuss the structure 
of lattice artifacts in the running coupling. To analyze which powers of the lattice 
spacing may occur, one has to list all terms which may appear in Symanzik's effective 
action p7| , ^5| , ^ 1 . In contrast to the standard situation (torus or thermodynamic limit), 
the presence of surfaces in the Schrodinger functional means that terms accompanied 
by one power of the lattice spacing are present in general. For our particular choice of 
boundary conditions 0, only the term 

-] 3 



a J d xtr F ok F ok \ 

x =0 

(3.2) 

and its image at xq = L are relevant. They can be cancelled by including one appro- 
priately chosen term in the lattice action, whose coefficient is denoted by Ct- When the 
proper dependence of ct on the bare coupling go is known, cutoff effects in the renor- 
malized coupling (and therefore also in -L max ) are reduced from 0(a) to 0(a 2 ). The 



perturbative expansion for ct(<?o) is known up to two- loop accuracy [3, |29| , |3C| ■ 

c t = 1 - 0.089c/g - 0.030c^ + . . . . (3.3) 

Considering that we have to insert values of g$ ~ 1, the perturbative series for this coeffi- 
cient appears rather useful. Nevertheless, as is always the case with perturbative expres- 
sions, it is an art to attribute an uncertainty to it (unless additional non-perturbative 
information is available). Below we will use both the full expression eq. (|3.3j ) and its 
truncation at one-loop order (i.e. ct = 1 — 0.089^) in our continuum extrapolation of 

-^max/^0- 



3.2.2 Continuum extrapolation of L max /ro 

The length scale, L mSuX , is defined as the value of L for which the running coupling, g(L) 
in the Schrodinger functional scheme has the specific value, 

5 2 (£max) = 3.48. (3.4) 

For a definite discretization this entails that L max /a is a unique function of (3, which 
can be determined by simulations of the Schrodinger functional Q . For ct to one-loop 
precision, numerical results are listed in Table 3 of that reference. We have repeated 
this calculation with ct given by eq. (|3.3| ) and increased precision, and our new results 
are presented in Table ||] below. 
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L/a f3 L/a (3 



4 


5.959(2) 


8 


6.476(3) 


5 


6.118(2) 


10 


6.654(3) 


6 


6.257(3) 


12 


6.799(3) 


7 


6.374(3) 


16 


7.026(4) 



Table 3: Bare couplings vs. lattice size at g (IS) = 3.48 and for two-loop c t as given in eq. (3.5). 



In order to form the ratio L max /ro we need both quantities L max /a and r$/a for 
the same values of (3. One of the data sets has to be interpolated. We used a linear 
interpolation (3 = Iq + l\ ln(L max /a), propagating the statistical errors. This linear 
interpolation is well justified: since the complete data sets can be fitted by adding 
a term, /2[ m (-^max/o)] 2 ) with a small coefficient, I2, a linear function in ln(L max /a) 
is locally an excellent approximation. As an alternative, we have also performed the 
mentioned global fit (with I2) and taken the fit as a representation of the data. The 
final conclusions remain unchanged. [] 

As regards the data for L max /Vo shown in Fig. [5], we note that the two-loop term 
—0.03 #o has quite a significant effect of the order of 8% at a/ro = 0.2, and its inclusion 
does indeed reduce lattice artifacts. On the other hand, some 15% of cutoff effects 
remain at a/ro = 1/5- We conclude that the perturbative expression eq. (|3.3|) is very 
useful, but due to the significant effect of the two-loop term we should not assume that 
the remainder can be neglected. 

We are led to model both data sets in a combined fit according to 



where ^ max 

ro 

tion 



U) a u) a? 

+ p) h p\ —!y , i= 1-loop, 2-loop , (3-5) 

ro ^o a =o r r£ 

and the 4 p's are fit parameters. This fit yields a continuum extrapola- 



-f'max 



ro 



= 0.718(16) . (3.6) 

a=0 



The figure shows nicely that the fit prefers a very small term linear in the lattice spacing 
after two-loop O(a) improvement. A value consistent with our final result eq. ( |3.6| ) is 
also obtained if one extrapolates just quadratically (the lower data points with two-loop 

3 When we take Ct to one-loop precision, the errors of L miH /ro are dominated by the errors in 
ff 2 (imax), while for our new data, the errors in ro/a are the larger ones. Therefore - despite the fact 
that the same values of ro /a enters both data sets - the data points shown in Fig. [5] are only weakly 
correlated. We neglect their correlation in the following. To be complete, we note that in fact a second 
source of statistical correlations of data points within each set is the interpolation (5 — lo + h ln(i max /a), 
where in most cases one simulation point contributes to two data points in Fig. |H| 
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Figure 5: Continuum extrapolation of L ma x/j"o- The two curves show the fit eq. { J3 - 5| ) who's value in 
the continuum limit is indicated by the square. 

O(a) improvement). However, the statistical error of such an extrapolation, which is 
five times smaller than the one quoted above, appears unrealisticaliy small. Only with 
a solid, non-perturbatively determined function ct(<?o) could we perform such a fit with 
confidence. 

4 Conclusions 

In this paper we have presented a precision computation of the low-energy reference scale 
tq in SU(3) Lattice Gauge Theory over a large range of f3, in which the lattice spacing 
varies by a factor of four. Our results at individual values of (5, which are displayed 
in Table ||, are supplemented by an interpolating parametrization, eq. ( |2.18| ), which 
provides estimates for r^/a at arbitrary values of (3 in the whole range of 5.7 < (3 < 6.57. 

Our study represents a significant improvement over previous calculations through 
the combined effects of variance reduction by means of the multi-hit technique, high 
statistics and the use of a variational approach. Finally, a precise definition of the force 
has been adhered to throughout our calculation. We observe a 2% deviation of our 
result at f3 = 5.7 compared to ref. §. 

Our results for the scale tq have been applied in the continuum extrapolation of 
the gap A of the static quark potential, and the low -energy scale L m ax) 

used to set the 

scale in the computation of the running coupling. In the case of A, the analysis revealed 
sizeable cutoff effects of 10-15% at lattice spacings of a ~ 0.1 fm. Nevertheless, a reliable 
continuum extrapolation could be performed, showing that at distances r ~ vq, the gap 
A is far from the universal value of bosonic string theories. 

The fact that ro/a has been calculated for small lattice spacings is crucial for a 
stable extrapolation of £ m axAoi since both o and a 2 lattice artifacts turn out to be 
relevant for this quantity. A reliable estimate for -L max /Vo is of importance in our study 
of the A-parameter, which will appear in a forthcoming publication |]26|| . 
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Our results for ro/a are presented in a form so that they can be easily used in 
any kind of scaling analysis of hadronic quantities. A typical example is the scaling 
behaviour of the vector mass, my tq for a fixed ratio of vector and pseudoscalar masses, 
e.g. my/mps = 0.7 [31, 3. Also, the value of m-proton^d hi the continuum limit can be 



obtained, which will serve to check how close ro is to the value 0.5 fm. 
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